Bulk and interfacial properties of binary hard-platelet fluids 
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Interfaces between demixed fluid phases of binary mixtures of hard platelets are investigated using 
density-functional theory. The corresponding excess free energy functional is calculated within a 
fundamental measure theory adapted to the Zwanzig model, in which the orientations of the particles 
of rectangular shape are restricted to three orthogonal orientations. Density and orientational order 
parameter profiles at interfaces between coexisting phases as well as the interfacial tension are 
determined. A density inversion, oscillatory density profiles, and a Fisher- Widom line have been 
found in a mixture of large thin and small thick platelets. The lowest interfacial tension corresponds 
to the mean bulk orientation of the platelets being parallel to the interface. For a mixture of large 
and small thin platelets, complete wetting of an isotropic-nematic interface by a second nematic 
phase is found. 
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PACS numbers: 61.20.Gy, 61.30.Hn, 68.05.-n, 82.70.Dd 



I. INTRODUCTION 



There is growing interest in thermal and structural 
properties of suspensions of platelike colloids like blood, 
clay sols, or liquid crystal dispersions because of possible 
applications in biomedicine (e.g. shape-selective separa- 
tion of cell components 0), geophysics (e.g. oil drilling 
0), or liquid crystal display technology [!|. From the- 
oretical considerations [j, 0| and computer simulations 
IE @ one expects an isotropic to nematic transition 
due to the orientational degrees of freedom for suffi- 
ciently high platelet concentrations. This transition is 
not observed for clays because the long-ranged Coulomb 
interaction between the charged clay particles leads to 
gelation Nonetheless, recent preparation meth- 

ods have been developed to produce suspensions of ster- 
ically stabilized platelets which do exhibit the expected 
isotropic to nematic transition 0, 0, fl3L IT3. It was 
shown experimentally and theoretically |l5l flrl ITtI 
IT^ | that polydispersity in the size of the platelets strongly 
affects the phase behavior. In particular, binary mix- 
tures of thin and thick platelets lead to an unexpected 
isotropic-nematic density inversion |l3lll6| . Although at- 
tention has been paid to the bulk phase behavior p^lis|. 
so far there are no studies on fluid-fluid interfacial proper- 
ties of such mixtures. On the basis of recent theoretical 
studies on fluids of thin hard platelets near hard walls 
|2C| , we expect that the competition of orientational 
entropy and excluded volume interaction leads to inter- 
esting fluid-fluid interfacial properties which should be 
experimentally accessible via optical techniques. 

In this paper, we study interfacial properties between 
coexisting isotropic and nematic phases of a binary mix- 
ture of hard square cuboids within the Zwanzig approx- 



imation |2l|. This model is chosen as the simplest non- 
trivial approximation to the above-mentioned properties 
of real platelet suspensions. Binary mixtures are used 
to mimic polydispersity. The hard-particle approxima- 
tion is doubtful for long-ranged platelet-platelet interac- 
tions, but there are model systems of quasi-hard platelets 
[Tllll^ ]. The actual shape of the platelets (disks, stripes, 
hexagons, or square cuboids) is expected to be less im- 
portant, only the thickness L and the width D are as- 
sumed to be relevant. Finally, the square face of the 
platelets can take only three rather than a continuous 
range of orientations (Zwanzig model |2l))- This model 
offers the advantage that the excess free energy functional 
can be determined within a fundamental measure theory 
|H |S1 and tnat the difficult determination of 
inhomogeneous density profiles becomes numerically fea- 
sible allowing one to study interfacial properties of binary 
hard-platelet fluids in detail. Due to the approximations, 
we only expect to find qualitatively correct results. 

This paper is organized as follows. In Sec. |n]we de- 
scribe the density functional and the fundamental mea- 
sure theory. Sec. IIIII presents two representative bulk 
phase diagrams of binary hard-platelet mixtures involv- 
ing isotropic, nematic, and columnar phases. In Sec. 
IIVI we determine the density and the orientational or- 
der parameter profiles as well as the interfacial tensions 
of isotropic-nematic interfaces in a mixture of large thin 
and small thick platelets. Sec. |3presents density profiles 
at nematic-nematic interfaces. The isotropic-nematic in- 
terfaces near an isotropic-nematic-nematic triple point 
are investigated in Sec. IVII Our results are summarized 
in Sec. IVTTI 
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II. DENSITY FUNCTIONAL AND 
FUNDAMENTAL MEASURE THEORY 

We consider a binary mixture of hard rectangular par- 
ticles of size Li x Di x flj, i £ {1,2}. The position of 
the center of mass r is a continuous variable, whereas 
the normal of the square face is restricted to directions 
(3 G {x, y, z}. The number density of the centers of mass 
of the platelets of size i and orientation (3 at position r is 
denoted by Qi.p{r). In the absence of external potentials, 
the equilibrium density profiles of the mixture minimize 
the grand potential functional 

V[{Qi,0}} = E / d3r ( kBT (Mft,/sWA|) - 1) 

F ex [{Qi,0}}, (1) 



i,0 



where the chemical potential fa and the thermal de 
Broglie wavelength Aj of platelets of size i are indepen- 
dent of the particle orientation. 

In fundamental measure theory, one postulates the fol- 
lowing form of the excess free energy functional [22l I23I 



F cx [{g t ,p}} = k B T / d 3 r$({n ; (r)}) 



with the reduced excess free energy density 



= -n ln(l-n 3 )+ J2 T 

cr£{x,y,z} 

(l-n 3 ) 2 
and the weighted densities 



n 3 



dV ^Jv-v')g^(v') 



(2) 



(3) 



(4) 



(5) 



for I e {0, lx, ly, Iz, 2x, 2y, 2z, 3}. 

The weight functions wj'I are obtained by expressing 
the Fourier transform of the Mayer-/ function as a sum 
of products of single particle functions. The Mayer-/ 
function equals —1 if the particles overlap and is zero 
otherwise. Explicit expressions of the weight functions 
are documented in Refs. [H El EU- In the limit 
of infinitely thin platelets (Lj/Dj — > 0) the excess free 
ener gy d ensity reduces to a third-order virial approxima- 
tion |20| . For convenience, we introduce the dimension- 
less variables fi* := fa/ksT + 31n(Di/Aj). 



III. BULK PHASE DIAGRAMS 

The bulk phase behavior of spatially homogeneous bi- 
nary platelet fluids has recently been studied within the 



Zwanzig model |17|. Rich phase diagrams are found in- 
volving an isotropic and one or two nematic phases char- 
acterized by different concentrations. The phase dia- 
grams are very sensitive to the size ratio D1/D2 and to 
the aspect ratios L\jD\ and L^jD^- 

Fig. ^shows two representative examples of bulk phase 
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FIG. 1: Bulk phase diagrams of binary mixtures of large 
and small hard platelets (L1/D1 = 0.01; (a): D1/D2 = 2, 
L2/D2 = 0.16; (b): D 1 /D 2 = 10, L 2 /D 2 = 0.1; see also the 
schematic side views Li x Di on the right side bar) in terms 
of the reduced chemical potentials /x* and /j, 2 . Full lines rep- 
resent strongly first-order transitions, whereas dashed lines 
denote second-order or weakly first-order transitions. In (a), 
one isotropic phase I, one nematic phase N, and one colum- 
nar phase C are found. The isotropic-nematic transitions in 
the monodisperse limits occur at fit = 0.829(1) for the large 
platelets (/zjj — > —00) and at fi 2 — 5.505(1) for the small 
platelets (/i* — » -co). Density and order parameter profiles 
for isotropic-nematic interfaces at state points a, b, and c are 
shown in Fig. |5] The nematic phase N demixes into two 
nematic phases Ni and N 2 . Fig. 0] exhibits density profiles 
for nematic-nematic interfaces at state points d and e, close 
to the lower critical point of the N\-N 2 demixing. For small 
values of nt, the nematic-columnar transition is of second or 
weakly first order (dashed line), whereas it is strongly first- 
order (solid line) if fj,* is sufficiently large. The dashed line 
appears to be straight because for this range of the chemi- 
cal potentials, the density of the large platelets is negligible 
compared to that of the small platelets. In (b), only a small 
part of the phase diagram around the triple point tri is shown, 
where the isotropic phase / coexists with two nematic phases 
ATi and N 2 . 

diagrams of binary mixtures of large and small platelets 
in terms of the reduced chemical potentials /z* and 
The coexistence curves in these diagrams are calculated 
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numerically by finding pairs of equilibrium states with 
equal chemical potentials and pressures. The phase tran- 
sitions at full curves are strongly first-order, whereas 
dashed lines denote second-order or weakly first-order 
transitions. 

For a binary fluid consisting of large thin and small 
thick platelets (see Fig. (Ha)), a single nematic phase 
TV is found, which is separated from the isotropic phase 
/ by a first-order phase boundary line. At larger val- 
ues of the two chemical potentials, the nematic phase 
N demixes discontinuously into two nematic phases N\ 
and N2- Furthermore, the nematic phase N undergoes 
a phase transition into a columnar phase C, where par- 
allel columns of small platelets are surrounded by single 
large platelets with orientations perpendicular to the col- 
umn axes. For small values of fi\, the nematic-columnar 
transition is of second or weakly first order. The phase 
boundary line appears straight because the density of the 
large platelets is negligible there, so that the mixture is 
nearly monodisperse. The nematic-columnar boundary 
line changes to strongly first-order at the tricritical point 
where the nematic-nematic boundary line ends. 

Upon decreasing the size of the small platelets, the 
lower critical point of this first-order N1-N2 coexistence 
curve shifts to smaller values of the chemical potentials 
until the N1-N2 and the I-N coexistence curves start to 
intersect, giving rise to a triple point at which two ne- 
matic phases N\ and N2 coexist with the isotropic phase 
I (see Fig. [T]»). 



IV. ISOTROPIC-NEMATIC INTERFACE 

In this section we discuss the properties of interfaces 
between the coexisting isotropic and nematic phases de- 
picted in the phase diagram of Fig. Qfa). The interface 
normal is always arranged to point in z-direction. The 
results are expressed in terms of the orientationally av- 
eraged number density 



Qi(z) := Qi, x (z) + Qi,y{z) + Qi, z (z), 
the nematic order parameter 

, x ._ QiA z ) - \{QiA z ) + Qj,y( z )) 
Si{z) :— . 

Qi{z) 

and the biaxial order parameter 

QiA z ) - QiA z ) 



Qi(z) 



Qi{z) 



(6) 



(7) 



(8) 



Fig. El shows the density and the nematic order pa- 
rameter profiles of the small thick platelets for the state 
points denoted as a, b, and c in Fig. ^a). The interface 
normal and the director of the nematic phase are parallel 
and the position z — of the origin is fixed by choosing 
Si(Q) — \- No biaxiality has been observed, that is, 
qi(z) = 0. We have confirmed that the density profiles 
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FIG. 2: Orientationally averaged density profiles Q2 (solid 
lines) and nematic order parameter profiles S2 (dashed lines) 
of the small thick platelets for the state points a, b, and c in 
Fig. 0^a) at a planar interface between the coexisting isotropic 
phase I at z — > 00 and the nematic phase N at z — * —00. The 
position z = of the origin is chosen to be given by S2(0) = 5. 
The interface normal and the director of the nematic phase 
are parallel, which corresponds to the lowest surface tension 
(see Fig. |3J. A density inversion occurs near the state point 
b. At the interface layers form which slightly prefer platelet 
orientations perpendicular to the interface normal. 



for both thin and thick platelets show oscillations of the 
same period, but the amplitude for the thin platelets is 
considerably smaller than for the thick platelets. 

Whereas the density of the thin platelets is always 
smaller in the isotropic phase than in the nematic phase 
along the coexistence curve, a density inversion for the 
thick platelets is found. For state b, the density of the 
small thick platelets in the isotropic and the nematic 
phase are nearly identical. This remarkable phenomenon 
of isotropic-nematic density inversion has been observed 
experimentally for a dispersion of sterically stabilized 
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gibbsite (A1(0H) 3 ) platelets 13]. 

For isotropic-anisotropic coexistence in binary mix- 
tures of long and short thin rods, one observes a large 
fraction of long rods in the anisotropic phase if there is 
a small fraction of the long rods in the isotropic phase 
[2^| . For platelets this fractionation does not occur to the 
same extend because the density g\ increases by at most 
one order of magnitude upon traversing the transition 
from the isotropic to the nematic phase. 

The nematic order parameter profiles shown in Fig. [21 
exhibit a small negative minimum (s2(z) < 0) on the 
isotropic side of the interface implying a depletion of 
platelets with orientation parallel to the interface. A sim- 
ilar minimum has been found for freely rotating platelets 
near a hard wall [l9j . This indicates that the nematic 
phase acts like a wall creating a corresponding depletion 
zone in the adjacent isotropic phase. 

The surface tension 77 at of an isotropic-nematic inter- 
face with the director of the nematic phase at z — > — 00 
parallel (||) or perpendicular (_!_) to the interface nor- 
mal is shown in Fig. |3| as a function of the pressure p 
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FIG. 3: Isotropic-nematic interfacial tension 7/ n of a binary 
mixture of large thin and small thick platelets (see Fig. ffl^a)) 
as a function of the pressure p along the line of phase coexis- 
tence. Through the choice of the boundary condition for the 
orientation of the nematic order at z — > —00 the director of 
the nematic phase can be parallel (||) or perpendicular (_L) to 
the interface normal. The monodisperse fluids, consisting of 
either the large or the small platelets, correspond to the end- 
points of the curves at the left and the right, respectively. The 
density profiles and the orientational order parameter profiles 
for the thermodynamic states marked a, b and c are shown in 
Fig. H 

along the line of coexisting phases. For all pressures, the 
parallel configuration has the lower surface tension by 
a factor of rather precisely three. As is apparent from 
Fig. |3 the surface tension is a non-monotonic function 
of the pressure and exhibits a maximum close to state c. 
The maximum value of 7/tv is about 20 times larger than 
the ones for the monodisperse fluids corresponding to the 
endpoints of the curves. This implies that the interface 
of the binary fluid is rather stiff as compared to the ones 
of the corresponding monodisperse fluids. 

Calculating the typical surface tension at the inversion 
point b (jj N Dl/kBT w 0.1) for the above mentioned 



sterically stabilized gibbsite platelets (Di w 170nm) at 

room temperature T = 300K yields yj N « 14nN/m, 
which perfectly agrees with the estimate provided in Rcf. 

0' . . 

A similar behavior has been found for binary mixtures 
of thin and thick rods 12711 . 



V. NEMATIC-NEMATIC INTERFACES 

We now turn our attention to the properties of inter- 
faces between two coexisting nematic phases whose cor- 
responding phase boundaries are depicted in the phase 
diagram of Fig. ^a). Fig. 0] shows the density pro- 
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FIG. 4: The density profiles gi (solid lines) and Q2 (dashed 
lines) of the large thin and small thick platelets of the bi- 
nary hard-platelet mixture, respectively, at the planar inter- 
face between two coexisting nematic phases (iVi and N2) for 
the states d and e indicated in Fig. 0^a). 

files of the large and small platelets for such interfaces. 
Pronounced oscillations are found for the state d (see 
Fig. on that side of the interface which is rich in 

the small thick platelets. The wave length of the os- 
cillations A = 0. 175(5)I?i is more than twice the thick- 
ness of the thick platelets. The exponential decay length 
£ = 0.131(4)Z?i is of the same order of magnitude. The 
wave length A and the exponential decay length £ are the 
same for the large and the small platelets. 

For the state e, closer to the lower critical point of 
nematic-nematic coexistence, there is no sign of oscilla- 
tions on either side of the interface. Upon approaching 
the lower critical point along the nematic-nematic co- 
existence curve, the density profiles Q\ and Qi turn into 
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monotonic functions of z and the interfacial width broad- 
ens. In simple fluids, the disappearance of oscillations in 
the density profiles is connected to the existence of the 
so-called Fisher- Widom line [2^, which divides the bulk 
phase diagram into regions where the asymptotic decay 
of bulk two-point correlation function is either monotonic 
or exponentially damped oscillatory [29l l3fj| ; this behav- 
ior of the two-point correlation function carries over to 
the density profiles. Recently, such Fisher- Widom lines 
have been found for a binary Gaussian core model [3j| 
and a model colloid-polymer mixture |32|. For the bi- 
nary platelet model under consideration, the determina- 
tion of the Fisher- Widom line requires the solution of 
Ornstein-Zernike equations for a six-component system, 
two species and three allowed orientations, which is be- 
yond the scope of the present study. However, from the 
observed vanishing of the oscillations of the density pro- 
files we conclude that there is a Fisher- Widom line and 
it intersects the coexistence curve between d and e. 

We have also calculated the surface tensions of 
nematic-nematic interfaces for different relative orienta- 
tions of the platelets and the interface. Again the small- 
est value is taken by the configuration with the mean bulk 
orientation of the platelets being parallel to the interface 
(compare Fig. |3J). Upon approaching the lower criti- 
cal point of nematic-nematic coexistence, the expected 
mean-field-like vanishing of the surface tension and of 
the N1-N2 density difference has been recovered. 



VI. ISOTROPIC-NEMATIC INTERFACE NEAR 
A TRIPLE POINT 
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FIG. 5: Equilibrium density profiles of large (a) and small (b) 
thin platelets (see Fig. flTb^l at various triple point undersatu- 
rationse := /u|-/x| tri £ {10 -3 , 5 x 1(T 4 , 10 -4 , 5 x 10 -5 , 1(T 5 } 
(from left to right) at I-N2 coexistence. The position z — 
of the origin is chosen as the midpoint of the bulk densities 
of Ni and N2 . The dashed lines represent the densities of the 
nematic phase Ni at the triple point. The isotropic phase I 
is at z — > 00, and the nematic phase N2 is at z — > —00. 



In this section we consider interfacial properties of 
coexisting isotropic and nematic phases near the triple 
point in the bulk phase diagram of Fig. ^b) . We choose 
e := [i\ — jLt| tn > at I-N2 coexistence as a mea- 
sure for the undersaturation with respect to the chem- 
ical potential of the small platelets at the triple point 
H\ tTi = 5.294376743251(1), n?" = 3.054579122021(1). 
Notice that for the present Zwanzig model the accuracy 
of the bulk coexistence data is higher than those one 
would be able to achieve for models of continuously ro- 
tating platelets because only a finite set of coexistence 
equations has to be solved instead of nonlinear integral 
equations. Hence, the undersaturation s can be defined 
rather precisely in the limit e -» within the present 
model. This is crucial for the subsequent wetting analy- 
sis and motivates the use of the Zwanzig approximation. 

Fig. [3] shows the density profiles of the large and the 
small platelets for various values of e. The density pro- 
files exhibit a plateau region, signaling the onset of the 
formation of the still metastable Ni phase. In the limit of 
e — > 0, the local density of the plateau region agrees with 
the one of the bulk phase Aq at the triple point, whereas 
the asymptotic densities at z — > ±00 are the ones of the 
coexisting phases / and N2- 

Fig. shows the interface thickness i as a function of 




FIG. 6: Thickness t of the I-N2 interface as a function of 
the undersaturation e := /12 — M2 trl from the triple point (see 
Figs, ^b) and|^Jl. The solid and the dashed line are cal- 
culated according to the inflection point and the midpoint 
method, respectively, as described in the text. The correla- 
tion length £ = 0.389(3)Di is determined from the slope of 
the logarithmic growth of t for e — > 0. 



the undersaturation e. There are various possible def- 
initions for the interface thickness. We have used the 
following two measures: 

• The inflection point method: The interface thick- 
ness is given by the distance between the outermost 
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inflection points of the density profiles. 



VII. SUMMARY 



• The midpoint method: The distance between the 
midpoints of the density profiles for the I-N± and 
the N1-N2 interface defines the interfacial thick- 
ness. 

Upon approaching the triple point, we find that t di- 
verges for e -> as t ~ -£lne + ii with £ = 0.392(4)£>i, 
t\ = — 0.92(4)Di for the inflection point method and 
£ = 0.386(2)Di, ti = -0.94(2)£>i for the midpoint 
method (see Fig. |SJ. The logarithmic divergence of t 
is consistent with complete wetting of the I-N2 interface 
by a Ni film at the triple point in the absence of alge- 
braically decaying interaction potentials 33]. Moreover, 
we have confirmed that the value of £ determined from 
the thickness of the interface is in agreement with the 
correlation length of the wetting phase N% at the triple 
point. The latter has been obtained from the asymp- 
totic exponential decay of the density profiles from I-Ni 
or N1-N2 interfaces towards their Ni bulk values at the 
triple point. 

The importance of interface fluctuations is determined 
by the so-called wetting parameter [34| 



u> 



47T7/AT 2 £ 2 



(9) 



which for the present system takes a value of order unity 
at the triple point. Fluctuations are therefore neither 
negligible nor dominant. 

In order to validate complete triple point wetting by an 
alternative method, the vanishing of the dihedral angle 
i9jVi of a lense of phase N\ at the triple point has been 
checked using the relation 



cos$ 



(7^,) 2 -(7^ 1 ) 2 -(7^ a ) ! 
iN 2 



(10) 



The interfacial tensions 7}^ = 0.05003056658fc B T/Df 
and 7$^ = 0.05600483075fc s T/£>2 are obtained by ap- 
proaching the triple point along the I-Ni and the N1-N2 
coexistence curves, respectively, from below in terms of 
the chemical potential of the small platelets [i\ (see Fig. 

Recently, complete wetting of an I-N2 interface by a 
Ni film has been predicted for a binary mixture of long 
thin and thick rods (L > D) with L\ = L 2 ,Di/D 2 > 4 
[35|,|3(|. The ratio £/L\ in those systems is of the same 
order of magnitude as the value £/D\ = 0.389(3) we have 
obtained for a binary mixture of large and small platelets 
(L < D) with Li =L 2 , D x /Ih = 10. 

At the emerging I-N\ interface, the orientational or- 
der parameter profiles s\ and s 2 exhibit slightly negative 
minima indicating the depletion of platelets with orien- 
tation parallel to the interface. The absolute values of 
these minima are two orders of magnitude smaller for s 2 
than for S\. 



We have studied bulk and interfacial properties of bi- 
nary hard-platelet fluids using density-functional theory. 
The platelets form square parallelepipeds with orienta- 
tions restricted to three mutually perpendicular direc- 
tions. A fundamental measure theory is used to derive 
the excess free energy functional. The grand potential 
functional is minimized numerically and phase diagrams, 
density profiles, orientational order parameter profiles, 
and surface tensions are determined leading to the fol- 
lowing main results: 

1. The bulk phase diagrams for two representative ex- 
amples of binary hard-platelet fluids (Fig. ^| in- 
volve an isotropic and one or two nematic phases of 
different concentration as well as a columnar phase. 

2. For a binary platelet mixture of large thin and small 
thick platelets, the density profiles of thick platelets 
at the isotropic-nematic interface exhibit oscilla- 
tions (Fig. [21). A density inversion for the thick 
platelets is found. Near the thermodynamic state 
marked as b in Fig. |2 the densities of the thick 
platelets in the isotropic and the nematic phase 
are nearly identical. The interfacial tension for 
platelets with their mean bulk orientation paral- 
lel to the interface is smaller than for the corre- 
sponding perpendicular configuration. The bidis- 
perse mixtures have significantly larger interfacial 
tensions than the corresponding monodisperse flu- 
ids (Fig. EJl . 

3. Both the large platelet and the small platelet den- 
sity profiles exhibit pronounced oscillations on one 
side of the nematic-nematic interface, provided the 
chemical potentials are sufficiently high (Fig. 21 
state d) . Upon approaching the lower critical point 
of nematic-nematic demixing the oscillations vanish 
and the interface broadens (Fig. 0] state e). 

4. Just above the isotropic-nematic-nematic triple 
point shown in Fig. IHb), the I-N2 interfacial pro- 
files thicken due to the formation of the incipient 
nematic phase N% at this interface (Fig. [5J. Fig. 
El demonstrates that the thickness of the wetting 
film diverges logarithmically upon approaching the 
triple point. Complete triple point wetting is con- 
firmed explicitly by observing the vanishing of the 
dihedral angle at the triple point. 
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